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For a noncentrosymmetric superconductor such as CePtsSi, we consider a Cooper pairing model 
with a two-component order parameter composed of spin-singlet and spin-triplet pairing compo- 
nents. We calculate the superfluid density tensor in the clean limit on the basis of the quasiclassical 
theory of superconductivity. We demonstrate that such a pairing model accounts for an experimen- 
tally observed feature of the temperature dependence of the London penetration depth in CePtsSi, 
i.e., line-node-gap behavior at low temperatures. 



o 

u 



X3 

o 



> 

o 

o 

-I— > 



I 

O 

o 



X 



I. INTRODUCTION 

Much attention has been focused on the superconduc- 
tivity in systems without inversion symmetry {e.g., Refs. 
[1I I2I3I I4 5 6, and references therein). Recently, CePtsSi 
was found to be a heavy fermion superconductor without 
inversion symmetry in the crystal structurei2iSi2ii£ This 
motivates more detailed studies of the superconductivity 
in noncentrosymmetric systems. The lack of an inver- 
sion center in the crystal lattice induces antisymmetric 
spin-orbit coupling^^'^^ responsible for a mixing of spin- 
singlet and spin-triplet Cooper pairings." In CePtsSi, this 
mixing of the pairing channels with different parity may 
result in unusual properties of experimentally observed 
quantities such as a very high upper critical field Hc2 
which exceeds the paramag netic limit)'^i^i^i^°i" and the 
simultaneous appearance of a coherence peak feature in 
the NMR relaxation rate T^^ and low-temperature pow- 
erlaw behavior suggesting line-nodes in the quasiparticle 
gg^p8,9ji0ji4j^ji6 (-ggg g^jgQ Ref. The presence of line- 
nodes in the gap of CePtaSi is also indicated by mea- 
surements of the thermal conductivitjJ^ and the London 
penetration depth. 

In CePtsSi, the superconductivity coexists with an 
antiferromagnetic phas o''i^i^i-^°i-^'^i-^^i-^^i^°i^^ (see also Ref. 
|2^ . Generally one may have to include this aspect when 
the low temperature thermodynamics is analyzed in this 
material. The London penetration depth, however, which 
is entirely connected with the superfluid density, contains 
exclusively the information on the superconductivity, and 
provides for this reason a very suitable probe of the low- 
energy spectrum of the quasiparticles associated with the 
superconducting gap topology. Experimental measure- 
ments of the London penetration depth o n polyc rystalline 
and powder samples are reported in Refs. llOligl We note 
that CePtsSi is an extreme type-II superconductor with 
the Ginzburg-Landau parameter k ~ 140,^'"'^° and the 
nonlocal effect can be safely neglected. 

For a noncentrosymmetric superconductor such as 
CePtaSi, we will consider a Cooper pairing model with an 
order parameter consisting of spin-singlet and spin-triplet 
pairing components. Based on the same pairing model. 



we previously investigated the nuclear spin-lattice relax- 
ation rate to explain peculiarities observed in y-^iSiSiiS^ 
In this paper, we calculate the superfluid density and 
demonstrate that this pairing model gives simultaneously 
an explanation of the powerlaw temperature dependence 
of the penetration depth at low temperatures in CePtaSi. 

This paper is organized as follows. In Sec. II, we de- 
scribe the electronic structure of the system without in- 
version symmetry and our pairing model. In Sec. Ill, the 
equations for calculating the superfluid density are for- 
mulated. The numerical results are shown in Sec. IV. The 
summary is given in Sec. V. In the appendix, we describe 
the derivation of the quasiclassical Green functions used 
to compute the superfluid density for the present pairing 
model. 



II. SYSTEM WITHOUT INVERSION 
SYMMETRY 

We base our analysis on a system considered in 
Ref. 's, where the lack of inversion symmetry is incorpo- 
rated through the antisymmetric Rashba-type spin-orbit 

COUphng3i'^6,27,28,29 



^9k ' ^Tii]' Cj^jj^kri' ■) 



(1) 



with 



(2) 



Here, & = (cTx, ^y,<yz) is the vector consisting of the Pauli 
matrices, (cfe,,) is the creation (annihilation) operator 
for the quasiparticle state with momentum k and spin rj. 
We use units in which Ti = = ^. a(>0) denotes the 
strength of the spin-orbit coupling. The antisymmetric 
vector Qk {g-k = —Qk) is determined by symmetry argu- 
ments and is normalized as (g|)o = li^^ k-p is the Fermi 
wave number and the brackets (• • • )o denote the average 
over the Fermi surface in the case of a = 0. 
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Generally we may classify the basic pairing states for 
a superconductor of given crystal symmetry, distinguish- 
ing the spin-singlet and spin-triplet statesi^ A general 
argument by Anderson''^ shows that the inversion sym- 
metry is a key element for the realization of spin-triplet 
pairing states. Hence, the lack of inversion symmetry as 
in CePtaSi may be detrimental for spin-triplet pairing 
states. In other words, the presence of the antisymmet- 
ric spin-orbit coupling would suppress spin-triplet pair- 
ing. However, it has been shown by Frigeri et al^ that 
the antisymmetric spin-orbit coupling is not destructive 
to the special spin-triplet state with the d vector paral- 
lel to gk {dk II gk)- Therefore, referring to gk given in 
Eq. we adopt the p-wave pairing state with paral- 
lel d vector, dk = A{—ky,kx,0)^ Here, the unit vector 
k = {kx,ky,kz) = {cos (j) sin 9, sin (j) sin 9, cos 9). A further 
effect of the antisymmetric spin-orbit coupling is the mix- 
ing of spin-singlet and spin-triplet pairing components^ 
Interestingly, only the s-wave spin-singlet pairing state 
(belonging to Aig representation of crystal point group) 
mixes with the above p-wave spin-triplet pairing state 
(for example, d-wave states cannot mix with this p-wave 
state because of symmetry) 

This parity-mixed pairing state is expressed by the or- 
der parameter. 



A(r,fc) 



*(r)o-o -hdfc(r) • & 



*(r)(To + A(r)(- kyCTx + kxO'y 



(3) 



with the spin-singlet s-wave component ^{r) and the d 
vector dk{r) — A{r){—ky, kx, 0). Here, the vector r indi- 
cates the real-space coordinates, and ctq is the unit matrix 
in the spin space. While this spin-triplet part alone has 
point nodes, the pairing state of Eq. can possess line 
nodes in a gap as a result of the combination with the 



s-wave component 



23.29.33 



In this paper, we choose the 



The Eilenberger equation which includes the spin-orbit 
coupling term is given ag^^i^'^i^^i^^i^^i^^ 



ivpis) ■ Vg + [iujnf^ - agt ■ S - A,g] = 0, (5) 



with 



T3 



& 
o-*'' 



CTo 
-aoj ' 



-ayCrcFy, 



Ok 



f Ok^o 
g-feO-Q 



f 9k^o 
y -gfco-o J ' 



9k = \/ -{~ky,kx,0), 



A 

-At 



(6) 



(7) 



(8) 



(9) 



(10) 



Here, v-p{s) is the Fermi velocity, the variable s indicates 
the position on the Fermi surfaces, and the commutator 
[a, b] — db — bd. The Eilenberger equation is supple- 
mented by the normalization condition^^*^ 
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(11) 



isotropic s-wave pairing as \E' for simplicity. 



where 1 is the 4x4 unit matrix. Because CePtaSi is a 
clean superconductor^iiS we neglect the impurity effect. 

To obtain an expression for the superfiuid den- 
sity, we follow the procedure developed by Choi and 
Muzikarii242ii4iSii& We consider a system in which a 
uniform supercurrent flows with the velocity Vs , and the 
gap function has the r dependence 



III. QUASICLASSICAL FORMULATION 

We will calculate the superfiuid density on the basis of 
the quasiclassical theory of superconductivityi^^i2^i2& Fol- 
lowing the spirit of the theory)i2& in this study we assume 
Iv^ljlAIja <C Ef {sf is the Fermi energy). We consider 
the quasiclassical Green function g which has the matrix 
elements in Nambu (particle-hole) space as 




(4) 



A(r,fc) = A'(fc)e' 



i2Mvs-r 



*(7o + A^-kyax 



and accordingly 



i2Mvs-r 



A(r) = Ae 



i2Mvs-r 



(12) 



(13) 



The bare electron mass is denoted by M. The matrix 
elements of the Green function (0J are expressed as 



g{r,k,iujn) = g'{k,iujn), 



(14a) 



where ojn — TTT{2n+l) is the Matsubara frequency (with 
the temperature T and the integer n). Throughout the 
paper, "hat" • denotes the 2x2 matrix in the spin space, 
and "check" • denotes the 4x4 matrix composed of the 
2x2 Nambu space and the 2x2 spin space. 



f{r,k,iujn) = f'{k,iu}n)e 



f{r,k,iujn) = f'{k,iuJn)e' 



i2Mvs r 



'i2Mvs-r 



(14b) 



(14c) 
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-/{r,k,iujri) = g'{k,iujn)- 



(14d) 



The Eilenberger equation |(SJ| is rewritten in a form with- 
out the r dependence as 



[z(a;„ + (Z)f3~agfc.5-A',5'] =0, (15) 



with 



q = iAIvp{s) ■ Vs 




A' 



A' 

-A't 



(16) 



(17) 



(18) 



We obtain the following Green functions from the 
Eilenberger equation and the normalization condition 
(see Appendix), 



g = ffiCTi + .9ii<Tii, 

with the matrices ai and an defined b}iiiS24i 

CTiJi = ^ ("^o ± Sk ■ , gk = {-ky, kx,0)- 
Here, k = {kx,ky,0) = (cos 0, sin 0, 0), and 



91 



/i 



* + A sin 9 



911 



/ii 



Br 



* - A sin ( 



Bi Bii 
The denominators Bi and Bu are given as 
Bi = ±v/(c^«+<z)2 + |* + Asin0|2, 

Bu = ±V(t^n + g)2 + |*- Asingp, 



(19a) 
(19b) 



(20) 



(21a) 



(21b) 



(22a) 



(22b) 



and the signs in front of the square root are determined 
by the conditions 



sgn(Re{5i}) = sgn(Re{a;„}) , 
sgn(Re{gii}) = sgn(Re{w„}). 



(23a) 
(23b) 



The Green functions labeled by the indices I and II 
belong to the two distinct Fermi surfaces which are split 
by the lifting of the spin-degeneracy due to the spin- 
orbit coupling. The density of states on those two Fermi 
surfaces is different from each other in general. We define 
the density of states (the Fermi velocity) as A^iji (i'i,ii) 
on the Fermi surfaces I and II. We also define a parameter 



5 {—1 < 5 < 1) which parameterizes the difference in the 
density of states. 



Ni - Nil 

2Nn ' 



(24) 



where 2No = Ni + Nu. 

The supercurrent J is composed of the regular part, 
—ing, of the Green function in Eq. Q.^^ J is expressed 

by ■ 

J = J dsNF{s)vF{s)tr[ao{-iTTg)] 



ttT 

i 



T r 

- [Ni{vigi) + Nii{viigii) 



(25) 



where N-p{s) is the density of states at the position s on 
the Fermi surfaces, "tr" means the trace in the spin space, 
and the brackets (• • • ) denote the average over each Fermi 
surface. In Eq. H25() . we have referred to Eqs. (|14a|l and 
(|19a|l . In order to calculate the superfluid density tensor 
Pij [Ji = PijVsj), we expand giji in Eq. (|21a|l up to first 
order in Vs (or q) , and substitute them into Eq. . The 
expression for is then obtained as^ 



MttTY^ 



Ni 




viiVijl^ + A sin 6*1' 
'\(c^2 + |* + Asin0|2) 
wmwiij l* - Asin^p 



3/2 



I* - AsinSp) 



3/2 



(26) 



Now, to compute the superfluid density tensor pij in 
Eq. H26|l . we need to assume a model of the Fermi sur- 
faces. For the shape of the Fermi surfaces, we adopt the 
spherical Fermi surface for simplicity, and thus viji — 
f^iji^^ = Will (cos sin 0, sin (/) sin 0, cos (?) a-nd the Fermi- 
surface average (• • • ) — (l/47r) J^^ d<j) dOsinO ■ ■ ■ 
This spherical approximation is justified for CePtaSi, 
since measurements of Hc2 give evidence for a nearly 
isotropic mass tensoriiiii We furthermore set the two 
Fermi velocities equal;^ 



vi = Vll 



(27) 



From Eq. H26() , the components of the superfluid density 
tensor are obtained as 

Pxx = {Mv^No)2TrT ^ 

de 1^' + Asin^Psin^ei 



1 

X - 
2 



Ci{S) 



2 (cj2 + |$ + Asin6i|2) 



3/2 



CiiiS) 



de |^'-Asin6iP 



sm 



2 (-^2 ^_ 1^ _ A sin 6*1' 



,3/2 



(28) 



4 



, ^6* I* + A sin 6'P sin 6*0082 I 



(w2 + |* + Asin6l|2 



3/2 



^6* I A sin 61 1 2 sin 6* 0082 6* 
2 (w2 + |4'-A8in0|2)^/' 



(29) 

Pxy = Pyx = 0, Pzx — Pxz — Pzy = Pyz = 0, and Pyy = 

Paja: . Here, the weighting factors Ci ji in the case of the 
model of Eq. H27(l are given as 

Ci((5) = l + <5, Cii{6)^l-5. (30) 

At zero temperature, 

p = Pxx{T = 0) = Pzz{T ^ 0) 



(31) 



where we have utihzed a formula for an arbitrary func- 
tion F, limT^o 27rT F{uJn) = J dujF{uj). 

The gap equations for the order parameters 5* and A 
are given by2^'"'* 

[K{f+)+6K{f-) 

\u)ri I <UJc 

+ A,„ (sin ) + S\rn (sin 6*/+) , 



A = 



where 



ttT Y [At (sin 6I/_ ) +(5 At (sin ( 

I 1^71 I <i^c 

+ A„(/+)+5A,„(/_)], 

/± = 7, > 



(32) 



(33) 



(34) 



and LOc is the cutoff energy. The coupling constants As 
and At result from the pairing interaction within each 
spin channel (s: singlet, t: triplet). Am appears as a scat- 
tering of Cooper pairs between the two channels, which 
is allowed in a system without inversion symmetry.2^ 

In the limit T — Tc (Tc is the superconducting critical 
temperature), the linearized gap equations allow us to 
determine As and At by 



As 
At 



1 / 2 ttJ 

A„i I - — h — 

w V3i^ 4 



1 + 



1 



- A„, V 



nS 
T 



0<ri<([^c/7rT-l)/2 



2 

2n+l' 



^ = -r 

when the parameters Am and v are given. 



(35) 
(36) 
(37) 

(38) 
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(b) 



Fermi Surface II 



FIG. 1: Schematic figures of the gap structure on the Fermi 
surfaces, (a) On the Fermi surface I, the gap is I^P -I- Asin&|. 
(b) On the Fermi surface II, the gap is jvp — Asin6'|. In these 
figures, it is assumed that both ^ and A are real and positive, 
and * < A. 



IV. NUMERICAL RESULTS 

In this section, we will show the numerically evaluated 
results for the superfluid densities pxx and p^z- To calcu- 
late them in Eqs. H28() and H29|l . we need the temperature 
dependence of the order parameters and A. We will 
use the order parameters obtained from the gap equa- 
tions [Eqs. and for the parameters Wc — lOOTc 
and Am = 0.2. This is a representative set of parameters. 
Different choices would not lead to qualitatively different 
results as long as the gap topology is not altered. We 
have calculated the superfluid densities also for Am = 0.1 
and obtained qualitatively same results. 

The Green functions in Eq. (|A.18|) are substituted into 
the gap equations. When solving the gap equations, A is 
fixed to be real without loss of generality, resulting in a 
real ^' as well>^ Referring to Eq. (|2rbl) . we notice that the 
superconducting gaps are |^ -I- A sin 6*1 and \^ — A sin 6*1 
on the Fermi surfaces I and II, respectively. Such a gap 
structure can lead to line nodes on either Fermi surface 
I or II (as shown in Fig. Pl.i^iSSiS When the signs of 
and A are reverse to (same as) each other, gap nodes 
appear on the Fermi surface I (II). The relative sign is 
controlled by the parameters Am, (5, and v in the present 
formulation. In this paper, we choose such parameters 
as * > and A > 0, so that gap nodes {\^ - Asin0| = 
0) appear on the Fermi surface II (see Fig. Under 
this circumstance, we can obtain stable order parameters 
^I^ and A when the difference in the density of states 5 
defined in Eq. is set as 5 > —0.2 for Am =0.1 and 
0.2. For (5 > —0.2, stable 4' and A are obtained when the 
singlet-to-triplet components ratio v defined in Eq. I|38|l 
is set as ^ 0.3 (Am = 0.1) and v > 0.5 (Am = 0.2). 

In Fig. 121 we show in a low-temperature region 
the reciprocal square root of the superfluid densities 
1 / \/ Pxx{T) and 1 / ^/p^jT), which correspond to the 
London penetration depth Al(T). We set here the pa- 
rameter ly — 0.6, for which the gap nodes are line nodes 
on the Fermi surface Ilm^ Indeed, the data exhibit the 
T-linear behavior at low temperatures, indicating the ex- 
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1 I ' ■ • 1 

0.05 0.1 0.15 0.2 



T/Tc 

FIG. 2: Plots of the reciprocal square root of the superfluid 
densities vs. temperature, 1/y/pxx (a) and 1/y/pzz (b), for 
several values of the difference in the density of states S. The 
singlet-to-triplet components ratio v is set as u = 0.6. For 
comparison, dotted lines indicate those for a point-node gap, 
i.e., pure p-wave gap |Asin6'| = 0). 



istence of line nodes. For comparison, we also plot in 
Fig. 13 the same quantities calculated for a point-node 
gap (dotted line), which is contrasting well with the 
linc-nodc-gap behavior. The present results explain the 
experimentally observed T-linear behavior of Xi,{T) in 
CePtsSiiSii^ We note here that CePtaSi is an extreme 
type-II superconductoriiS and nonlocal effects can be 
neglected. 

Concerning the dependence on 6 the difference in the 
density of states [Eq. (|24ll ] , we notice in Fig. |21 that 
the smaller S {—1 < S < I), the stronger temperature- 
dependence appears. Note the nearly temperature- 
independent behavior of the state with point-nodes (dot- 
ted line) at low temperatures corresponding to a be- 
havior. This indicates weaker contributions of low-energy 
quasiparticles. For the identical Fermi velocities on the 
two Fermi surfaces, the smaller 5 {—1 < S < 1) leads 
to a smaller weighting factor Ci — 1 + 6 and a larger 
Cii = 1 — 6. Therefore, with decreasing 6, to the super- 
fluid densities the contribution of the fully-gapped Fermi 
surface I shrinks with the decreasing weighting factor 
Ci = 1 + 6. On the other hand, the contribution of the 
Fermi surface II with the gap nodes increases because of 




T/Tc 

FIG. 3: The temperature dependence of the superfluid den- 
sities, pxx (a) and pzz (b), for several values of the difference 
in the density of states S. The singlet-to-triplet components 
ratio f is set as = 0.6. 



the growing weighting factor Cn = 1 — 6 and the effect 
of the line-node-gapped Fermi surface II is enhanced by 
decreasing 5 (Fig.[2Jl. 

In Figs. 13 and ^ we show the superfluid densities 
Pxx{T) and Pzz{T) to see the overall temperature de- 
pendence. The quantities correspond to 1/\^{T). With 
decreasing 6 {—1 < (5 < 1) in Fig. [SJ the curves deviate 
from an upper convex curve and become gradually upper 
concave curves, namely they deviate gradually from fully- 
gapped s-wave behavior (i.e., pii{T = 0) — pii{T) ~ T'* in 
the s-wave case) because of the same reason mentioned 
above for the 6 dependence in Fig. |21 

We show in Fig. 0] the dependence on the singlet-to- 
triplet components ratio v [Eq. H38() ] . With increasing v 
in Fig. 2Ia), the curvature of the upper concave curves 
of pxx becomes larger. On the other hand, the v depen- 
dence of pzz is weaker than that of pxx^ as seen in Fig. 
0Jb). The quantity p^x [Eq. (|28|l ] senses the gap topol- 
ogy emphatically near the equator of the Fermi surfaces, 
while pzz [Eq. H29|l ] senses it near the poles. For the 
singlet-to-triplet components ratio v 0.5, the places 
(or the angles 9) of gap nodes at which — Asin^j = 
(see Fig. are sufficiently away from the poles, and ap- 
proach gradually to the equator on the Fermi surface II 
with increasing v. Therefore, p^x (fzz) is sensitive (not 
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FIG. 4: The temperature dependence of the superfluid den- 
sities, pxx (a) and pzz (b), for several values of the singlet-to- 
triplet components ratio v. The difference in the density of 
states is set as 5 = —0.2. 



sensitive) to the change of for > 0.5 as seen in Fig.^] 
It is interesting to note the difference in the tempera- 
ture dependence between p^x and pzz- In Fig. we plot 
the ratios Pzz/Pxx as functions of the temperature for 
several values of v. They exhibit a nonmonotonic tem- 
perature dependence in contrast with a monotonic one in 
the cases of the axial state (point-node gap at the poles) 
and the polar state (line-node gap at the equator) shown 
in Ref. B 

We have so far shown the results obtained for the 
Fermi velocity model of Eq. (|?7jl {vi — vu) and Eq. 
(|30|l (Ci,ii((5) = 1 ± (5). Alternatively, we have cal- 
culated the same quantities also for a different model: 
oc 1/Niji [thus, Vu = vi{l + 6) / {1 — S) instead of Eq. 
(|27|l ] . The weighting factors are accordingly Ci = 1 + (5 
and Cu = (1 + <5)^/(l — S) in this model, and p = 
Pxx{T - 0) - pzz{T = 0) = {2MvfNo/3)il + 5)/{l - 6). 
Because of these weighting factors Ci = 1 + S and 
Cii = {1 + 5)V(1 - S), the larger S (-1 < S < 1), the 
more enhanced the contribution of the Fermi surface II 
becomes. It is an opposite dependence on S as compared 
to the case of Eq. ^ {Ci,n{S) = 1±S). The results for 
the superfluid densities are qualitatively similar to those 
obtained for the previous model of Eq. H27|) (uj — vu) 
with replacing S ^ —S in Figs. |2HS1 We plot, never- 
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FIG. 5: Plots of the ratio pzz/pxx vs. temperature for sev- 
eral values of the singlet-to-triplet components ratio i/. The 
difference in the density of states is set as S = —0.2. 
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FIG. 6: The temperature dependence of the superfluid densi- 
ties pxx and Pzz for the Fermi velocity model, uiji cx 1/M,n. 
u = 0.6 and S = 0.9. 



theless, a result for an extreme case in Fig. El where we 
set 1/ — 0.6 and 6 — 0.9. In this case, the superfluid 
densities are predominantly determined by the contribu- 
tion of the Fermi surface II with gap nodes, owing to the 
extreme value S = 0.9 and the resulting weighting factor 
Cii = (1 + ^)V(1-^) > Ci = 1 + S. It is noticed in Fig.EI 
that the superfluid densities are suppressed at high tem- 
peratures. The temperature dependence of pxx in the 
region T < 0.3Tc of Fig. El is well fitted into an experi- 
mental result for 1/X'^{T) in CePtsSiji^ and an unusually 
strong suppression of l/Xf^ (T) at high temperaturesiS. is 
somewhat similar to that of pzz in Fig. El However, the 
difference in the density of states estimated from an band 
calculation for CePtsSi is \S\ - 0.25-0.3.4 Therefore, the 
strong suppression of 1/Xj (T) at high temperatures ob- 
served in CePtgSi (Ref. [13) remains to be accounted for 
at this moment. 



7 



SUMMARY 



APPENDIX 



We calculated the temperature dependence of the su- 
perfluid densities pxx{T) and Pzz{T) for the noncen- 
trosymmetric superconductor with the Rashba-type spin- 
orbit coupling represented by Eq. We showed that 
the gap function of Eq. Q , which has the spin-singlet and 
spin-triplet pairing components, explains the line-node- 
gap temperature dependence of the experimentally ob- 
served Al(T) in CePtaSiiiSiiS While the low-temperature 
behavior (T < 0.2Tc) of 1/A£(r) can be reproduced qual- 
itatively by that gap function, the high-temperature one 
still remains to be accounted for. 

The detailed information on the Fermi surfaces in the 
actual material CePtaSi is not available so fari^ The 
main difficulty here lies in the fact that this material is 
a heavy fermion system with strongly renormalized car- 
riers. If the Fermi velocity on the Fermi surface with 
gap nodes (on the Fermi surface II in our assumption) is 
sufficiently large in the case of the model wiji oc 1/Niji 
(otherwise, if Nu ^ Ni in the case of vi ~ vu), the 
anomalously strong suppression of l/Xl{T) observed ex- 
perimentally at high temperatures^^ may be explained as 
in Fig. El (where Cu ^ Ci). Such an assumption seems 
not unreasonable in view to the large renormalization 
factors for the effective mass suggested from thermody- 
namic measurements. On the other hand, the London 
penetration depth was not measured on a single crystal, 
but on polycrystalline and powder samples, ^"■"'^^ to which 
the anomalous behavior at high temperatures could be 
attributed. We also note that the superconducting tran- 
sition around Tc is rather broad in CePtaSi at least at 
this moment Thus also an unusual behavior of the su- 
perconducting phase close to Tc may play a role in the T- 
dependence of 1/X^{T). In any case, further experimen- 
tal studies (e.g., experimental measurements on a single 
crystal) and theoretical studies using more information 
on the Fermi surfaces involved in superconductivity are 
needed in the future to accomplish a detailed fitting. It 
would be also important to test experimentally CePtaSi 
for the intriguing nonmonotonic temperature dependence 
of Pzz/Pxx (see Fig. which could provide information 
on effective parameters of the model. 
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In this Appendix, we describe the procedure for de- 
riving the quasiclassical Green functions for a noncen- 
trosymmetric superconductor with the Rashba-type spin- 
orbit coupling represented by Eqs. Q and (0) and with 
the gap function of Eq. J^Jl. The explicit form of the 
Eilenberger equations which will be given here [Eqs. 
(jA.lOl) and (|A.13p would be useful for future studies in 
inhomogeneous systems such as surfaces, junctions, and 
vorticesSIS in the noncentrosymmetric superconductor. 

We start with the Eilenberger equation given in Eq. 
(0), namely^ 



idg + [iunfs - agu ■ S - K,g\ = 0, 



(A.l) 



where we have defined d = v-p{s) ■ V. It is convenient to 
define, for the Green functions [Eq. , 

(A.2) 

We also define, for the gap function [Eq. , 

A = koiay, At = -idykl, (A.3) 

where 

Ao(r, k) = *(r)ao + A(r) {-kyGx + hdy) . (A.4) 
The Eilenberger equation (|A.1|I is then written down to 
dgo + iagk ■ (o-go - goa) 

+ (Ao/o - /oAj) = 0, (A.5a) 

(5 + 2cj„)/o + iagk ■ (o-fo - fao-) 

+ (Aoffo - 5oAo) = 0, (A.5b) 

{d - 2w„)/o + iagk ■ (a- Jo ~ foO-) 

+ {Also - loAJ) = 0, (A.5c) 

dgo + iagk ■ - 50*^) 

-f (At/o - /oAo) - 0, (A.5d) 

where gk = ^/3/2{-ky, kx,0)- 

From Eq. the normalization conditions are 



9o + /o/o = ^0, 



9o + /o/o = o'o, 



5o/o — -/o5o, 



(A.6a) 



(A.6b) 



(A.6c) 
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ffo/o - ^/o.9o- 



(A.6d) 



Now, we consider a rotation in the spin space which is 
represented by the matrices {Ul.Uk — Uk&l = (^Q)r^^ 



~ 1 n -k 



u=^ 



1 k' 



V2\-k'- 1 



,(A.7) 



where k = (kx,ky,0) = (cos (/), sin 0, 0) and k'j. = ky ± 
ikx- A physical meaning of this rotation is that after the 
rotation the quantization axis in the spin space becomes 
paraUel to the vector gk \\ {—ky,kx, 0) at each position on 
the Fermi surface. We define the matrix elements 
of the Green functions after the rotation as 



ulg^Uk : 



ga gb 
gc gd 

la fb 
fc Id 



ulfoUk = 
ulloUk 



fa fb 
fc fd 

ga gb 
gc gd 



(A., 



We apply ul and Uk to the Eilenberger equations l|A.5|l 
from left and right respectively, so that we obtain the 
following sets of equations in the case when gk and Aq are 
given in Eqs. (O and l|A.4(l . Here, we define a' = 
and 

Ai = * + A sin 6*, An = * - A sin 6. (A.9) 
For the Green functions with the suffix "a", 

dga + Aufa - Ai*i/, = 0, (A.lOa) 

(5 + 2w„)/, + A„g, - A„g, = 0, (A.lOb) 

(d - 2uOn)fa + A*ig„ - A*ig„ = 0, (A.lOc) 

5ga + A*i/, - A„/, = 0. (A.lOd) 

For the "6", 

dgb - 2ia' sin%b + An fb - A^fb = 0, (A. 11a) 

{d + 2LJn)fb - 2ia' sine fb + Angb - Aigb - Q[A.llb) 

[d - 2un)fb - 2ia' sine fb + A^.gb - Algb = O^A.Uc) 

dgb - 2ia' sin dgb + Ai*i/b - Ajb = 0. (A. 1 Id) 
For the "c", 

dgc + 2iol sin dg, + Ai/e - A\^f, = 0, (A.12a) 
(d + 2w„)/c + 2iol SYndf^ + Ai5c - Aii.gc = C(A-12b) 
{d ~ 2lj„)/, + 2iol sine fc + A*g, - Ay^ = 0(,A.12c) 



dgc + 2m' sin % + A*fc ~ An fc = 0. (A.12d) 
For the "d", 

dgd + Aifd-Alfd^O, (A.13a) 

{d + 2iUn)fd + Aigd~Aigd = 0, (A.13b) 

{d - 2L^„)/d + Alga ~ A^ga - 0, (A.13c) 



dgd + A*Jd~Aifd = 0. 



(A.lSd) 



Note above that the Green functions with the suffixes 
"a", "b", "c", and "d" are decoupled from each other in 
these sets of the Eilenberger equations. 

We can solve Eqs. (|A.11|I and (jA.12p in the case of spa- 
tially uniform system, and find that gb = = = 5b = 
and gc — fc = fc — gc — Q for a' ^ 0. We also notice 
that the Green functions with the suffixes "b" and "c" 
are zero everywhere (even if the order parameters Ajji 
are spatially inhomogeneous), when solving the differen- 
tial equations IjA.lljl and (|A.12(I with the initial values 
equal to zero. On the other hand, the Green functions 
with the suffixes "a" and "d" in Eqs. HA.10|I and HA.13|I 
have finite values in general. Therefore, we can rewrite 
Eq. lEM as 



ulgoUk 
ulfoUk 
We obtain accordingly 



J' ^If^^^ 



fa 
fd 



ulkUk 



fa 
fd 

ga 
gd 



go = 



fo 



1 f _gd + ga 

2 \-kL{gd-ga) 



-k'+i9d - ga) 

gd + ga 



_fd + fa ~k'^ifd-fa) 
~k'_{fd-fa) fd + fa 



(A. 14) 



(A.15a) 



(A.15b) 



f fd + fa -k'+{fd-fa)\ .A1^p^ 

2\~k'_{fd~ fa) fd + fa ^^■^^^> 



9o 



gd + ga 



2 \-ki{gd - ga) gd 



~k'+i9d--ga)\ (A.15d) 



and finally, from Eq. (|A.2p . 



g = gi<Ji + giimi, 

f = {fl^l + /iio-ii)i(Ty, 

/ = -«0-,y (/l (71 + /llO-Il), 

9 = ^^y{gi^i + gii^ii)S-y 



(A.16a) 
(A. 16b) 

(A.16c) 
(A.16d) 
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with 



(To±gfe-(T , gk ^ {-ky,kx,0). (A. 17) 



Here, we have defined Ai = Ad and An = Aa, {A = 

{9JJ,g})- 

In the case of spatially uniform system, we obtain from 
the Eilenberger equations ljA.10|) and l|A.13|) and the nor- 
malization conditions in Eq. ||A.6|) . 



91 = gd = 



911 = 9a (A. 18a) 

mi 



h = fd — 



Ai 



/ll = fa = 



Bn 



(A. 18b) 



/l = fd 



91 = 9d 



Aj* - _ - 

/ll = fa 

Bl 



9ll = 9a 



with 



-Bi II — \l u, 



^11 
Bn' 



Br 



lAi.iil 



(A.18c) 



(A.18d) 



(A.19) 



At last, in these equations we replace — > w„ + g com- 
paring Eq. (|15|) with Eq. (|A.1|I . so that we obtain the 
Green functions in Eqs. ()21|l . 
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